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Abstract 

Multifragnientation (MF) results from lA GeV Au on C have been compared 
with the Copenhagen statistical multifragnientation model (SMM). The com- 
plete charge, mass and momentum reconstruction of the Au projectile was 
used to identify high momentum ejectiles leaving an excited remnant of mass 
A, charge Z and excitation energy E* which subsequently multifragments. 
Measurement of the magnitude and multiplicity (energy) dependence of the 
initial free volume and the breakup volume determines the variable volume 
parameterization of SMM. Very good agreement is obtained using SMM with 
the standard values of the SMM parameters. A large number of observables, 
including the fragment charge yield distributions, fragment multiplicity dis- 
tributions, caloric curve, critical exponents, and the critical scaling function 
are explored in this comparison. The two stage structure of SMM is used to 
determine the effect of cooling of the primary hot fragments. Average frag- 
ment yields with Z > 3 are essentially unaffected when the excitation energy 
is < 7 MeV/nucleon. SMM studies suggest that the experimental critical 
exponents are largely unaffected by cooling and event mixing. The nature of 
the phase transition in SMM is studied as a function of the remnant mass and 
charge using the micro canonical equation of state. For light remnants A < 
100, backbending is observed indicating negative specific heat, while for A > 
170 the effective latent heat approaches zero. Thus for heavier systems this 
transition can be identified as a continuous thermal phase transition where a 
large nucleus breaks up into a number of smaller nuclei with only a minimal 
release of constituent nucleons. Z < 2 particles are primarily emitted in the 
initial collision and after MF in the fragment deexcitation process. 
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I. INTRODUCTION 

Multifragmentation (MF) is the dominant decay mode in heavy ion reactions when the 
excitation energy is comparable to the nuclear binding energy. The recent use of Att detectors 
capable of observing a substantial fraction of the particles and fragments emitted in a given 
interaction has been invaluable in furthering the elucidation of this complex phenomenonETQ. 
Several review articles of MF have been published in recent yearsQu. 

One of the most complete MF experiments to date has been that performed by the EOS 
collaboration for the ^^"^Au on C systemll3nl3. The use of seamless detectors such as the EOS 
time projection chamber coupled with the multiple sampling ionization chamber (MUSIC 
II) permitted the observation of practically all the charged particles and fragments emitted 
in each event, ranging from protons to heavy fragments. Full reconstruction was therefore 
possible for a large fraction of the events. The high energy asymmetric collision of a lA GeV 
projectile on a light target is uniquely favorable for the kinematic separation of the initial 
nuclear collision and the subsequent MF transition. This data set permitted us to establish 
that MF occurs following the expansion of a remnant formed with charge Z, mass A and 
excitation energy E* after the prompt ejection of high momentum light particlesO'Ej. We 
have applied model independent methods used in the study of critical phenomena to extract 
the value of several critical exponentsE^tJ. The critical scaling function, which describes the 
behavior of the nuclear remnant near a critical point, was determinedllj. In performing this 
analysis we were aided by percolation calculationstatJ, which served as a quantitative guide 
to the application of the methods developed for the study of critical exponents in small 
systems. 

The EOS experiment was an outgrowth of earlier inclusive studies of MF in the inter- 
action of Xe and Kr with high-energy protonscHca. In this work high precision counter 
techniques were used to obtain accurate information about the kinetic energy spectra of 
isotopically resolved nuclear fragments. The reduced Coulomb barrier seen in the fragment 
spectra indicated that nuclear fragments are emitted from an expanded nuclear system. Sys- 
tematics of the fragment kinetic energy spectra also showed that the remnant was lighter 
than the target nucleus implying that the nuclear remnant existed for a time after the initial 
collision. The relative fragment yields of 63 isotopes could be understood using a thermal 
droplet modelE3 with a free energy parameterization based on the semi-empirical mass for- 
mula. For energy depositions of ~8 MeV per nucleon this multi-isotope thermometer gave 
a freezeout temperature of ~5 MeV. The systematic variation of the "kinetic temperature" 
(exponential slope parameter extracted from fragment kinetic energy spectra) as a function 
of the fragment mass also indicated that Fermi momentum in the remnant must play a 
significant role in the fragment formation process. 

More speculative was the suggestion that the power law yield of the fragment masses 
Y=A~'^ with T ~2.6 observed in 50-400 GeV p on Xe collisions showed that nuclear frag- 
mentation might give information about a possible "liquid-gas" thermal phase transition in 
charged nuclear matterBil. Important support for this view came from the exclusive emul- 
sion data of Waddington and FreieiS, which were analyzed by Campi to show that the 
conditional moments of the individual fragment events exhibited features characteristic of a 
critical transitioned]^. The occurrence of thermal equilibrium in the MF process was strongly 
supported by the experiments of the ALADIN collaboration, which showed that fragment 



yields were independent of the entrance channel when the data was scaled for projectile or 
target niassoEl. 

The above cited evidence for the formation of a substantially equilibrated remnant which 
expands prior to MF suggests that it would be appropriate to compare the EOS results with 
a thermal model in which MF occurs from an expanded state. Two widely used thermal 
multifragmentation models are the Copenhagen (SMM)eHI and Berlin (MMMC)E3 statistical 
treatments. These models differ in their parameterization of the expansion and in technical 
details but are similar in their thermodynamic approach. For very small systems both 
these micro canonical models predict the onset of a very inhomogeneous state at a definite 
excitation energy. This inhomogeneous excited state consists of a number of normal density 
fragments accompanied by a statistically insignificant number of nucleons. Subsequently 
these fragments cool by light particle emission. We present here a comparison of the EOS 
data with SMM. Several data and model comparisons have been made previouslyE3c3. These 
comparisons seriously suffered from the fact that the experiments did not determine the Z, 
A, and E* of the remnant. Instead, these Z, A, E* values, which constitute the proper input 
to test models of MF, were obtained either by use of a dynamical first stage modelEj'E^, by 
using a comparison of some of the data with SMM to constrain theniSED&cll, or by more 
complex backtracing from fragment datafi^cHl. A great advantage afforded by asymmetric 
reverse kinematics collisions is that they permit an accurate separation of the initial reaction 
phase from the subsequent decay of the excited remnant. 

The use of SMM also permits us to examine several aspects of the EOS critical exponent 
analysis that are potentially problematic. The EOS results were obtained for fragments in 
their final, cold state. However, the fragments are presumably formed in an expanded hot 
state. In SMM the remnant is equilibrated, the fragments are formed in the hot system, 
then separate under the influence of the Coulomb force and undergo deexcitation. As has 
been noted elsewherec^, the distribution of the hot, primary fragments may be affected by 
secondary decay. The difference between the two distributions could, in principle, affect 
the values of the critical exponents. Since SMM gives separate results for hot and cold 
fragments, the effects of secondary decay can be probed. 

The EOS results indicate that the first prompt step leads to a distribution of remnants. 
The analysis groups these remnants according to the total charged particle multiplicity, 
m, which serves as the control parameter. Events characterized by a given multiplicity 
will generally include a range of remnant Z, A, E* values. It has been noted that such 
event mixing can affect the values of the critical exponentsea. Furthermore, the use of 
multiplicitv instead of temperature or excitation energy as the control parameter has been 
quest ionedE3c3. Comparison of SMM results obtained for the experimental distribution of 
remnants with those obtained for a single remnant permits us to probe the importance of 
these effects. 

This paper is organized in the following manner. Section II gives a brief summary of SMM 
and shows how the data determine the variable volume of SMM. Only a single adjustable 
parameter remains. Section III summarizes the properties of the experimental remnant 
distribution, which serves as the input data to SMM. The occurrence of radial expansion 
energy as part of the excitation energy is discussed in this section. The comparison between 
experimental and model results is presented in Section IV. The various factors that can 
affect the fragment yields, the extraction of critical exponents, and the critical multiplicity 



are considered in Section V. In Sec. VI, the physics of MF and the nature of the thermal 
phase transition in SMM as a function of the remnant mass and charge is explored. The 
identification of the order of a thermal phase transition for the MF of Au on C requires the 
use of the microcanonical equation of state. A summary of the results and our conclusions 
are given in Section VII. 

II. THE SMM MODEL 

SMM is a statistical description of the simultaneous breakup of an expanded excited 
nucleus into nucleons and hot fragmentsE^. Individual fragments at normal nuclear density 
are described with a charged liquid drop parameterization. The free energy of a fragment 
A, Z {Z > 3) is given by 

^ A,Z ^ tr-ans ~r i^ vol ~r ^ surf ~r ^ sym ~r ^ Coul \^) 

and is used to determine the fragment formation probability. This solution explicitly as- 
sumes the inhomogeneous nature of the hot MF final state. Light fragments Z < 3 may 
also be present in the hot MF final state. For the Z >?> fragments, a quantum mechanical 
description is used for the temperature dependent volume, surface, and translational free 
energy of the fragments. The temperature independent parameters are based on the coef- 
ficients of the semiempirical mass formula. The critical temperature, at which the surface 
tension of neutral nuclear matter droplets would go to zero, is in the range suggested by 
infinite neutral nuclear matter calculationseS. 

In SMM the translational free energy depends on the free volume. The free volume, Vf, 
can be expressed in terms of the volume of the multifragmenting system at normal nuclear 
density, Vrem, 

Vf = xVrem (2) 

where the free volume parameter x depends on the SMM fragment multiplicity according 
to the relation: 
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where Ro=l-l7 A^^^ fm and M is the charged plus neutral hot fragment multiplicity. The 
crack width parameter, d, scales the magnitude of the multiplicity dependent free volume. 
The breakup volume VJ,, which includes the volume of the fragments, is VJ, = (1 + K)Vrem, 
where n is the Coulomb reduction parameter^Hl. 

We have previously shown that energy deposition in the Au on C reaction is proportional 
to nucleons knocked out of the Au nucleus by quasi-elastic nucleon-nucleon scatteringli3. 
If we assume that the excited remnant initially is produced in the Au volume, then the 
experimentally determined initial free volume is given by: 

V} = Vau (Aau - Arem)/AAu (4) 

and the remnant volume is given by: 

6 



K-em — yAu [^rem/^Au) (5j 

Proton-proton correlation experiments for the multifragmentation of lA GeV Au on Au 
show that the freeze out volume is ~2 V^^, nearly independent of excitation energy and 
remnant mass and charged. 
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FIG. 1. a) SMM free volume (Eq. 2) and the experimental initial free volume (Eq. 4) for the 
MF region, b) SMM breakup volume and the experimental freeze out volume from p-p correlations. 
The variation in V^ is due to the decrease of the remnant mass. The remnant does not expand in 
compound nucleus events which occur at low multiplicity. 

Fig. la shows the SMM free volume Vj, and the initial free volume V} as a function of the 
multiplicity. As expected, the SMM free volume Vf is about twice the initial free volume, 
consistent with the expansion from Vau to 2Vau- Note that the slopes of the Vf and V^ 
versus multiplicity curves track over the MF region. 

Fig. lb compares the breakup volume 

Vb = (1 + l^)Vrem = (1 + l^){Arem/ ^Au) Vau (6) 

with the freeze out volume from p-p correlations as a function of the multiplicity. The 
breakup volume Vf, tracks the freezeout volume from the p-p correlation experiments. The 
freeze out volume calculated in our earlier publicationll2l is different as it was obtained using 
initial volume of Au nucleus and not the remnant volume. Thus, experiment confirms the 
structure (M dependence) and scale (crack width parameter d) of the volume parameter- 
ization of SMM for the lA GeV Au on C experiment. In addition to these experimental 
arguments, theoretical arguments based on BUU modeling of the nuclear collisions also 
suggest that the expansion of the remnant is energy dependent. 

The remaining parameter is eo, the inverse level density parameter. The eo values were 
determined by comparison of SMM with the various experimental fragment yield distribu- 
tions and it was found that eo = 16 MeV. Thus, the so-called standard values of all these 
parameters turn out to give the best agreement with the data. (See Table I.) 

In SMM the primary fragments are propagated in their mutual Coulomb field and then 
undergo secondary decay. In the model successive particle emission from hot fragments 



with A>16 is assumed the deexcitation mechanism. The deexcitation of these fragments 
is treated by means of the standard Weisskopf evaporation model. Light fragments (A < 
16) deexcite via Fermi breakup. The hghtest particles (A<4) can be formed only in their 
ground states and undergo no secondary decay. We have used a version of the model that 
incorporates only thermal degrees of freedom. Consequently, radial expansion or angular 
momentum are not included in this comparison between data and SMM. 

III. THE INPUT DATA 
A. Properties of the remnant 

The reverse kinematic EOS experiment permitted the identification of charged projectile 
fragments 1 < Z < 79 in lA GeV ^^^Au on C interactions with high efficiency. The momenta 
of these fragments were measured and used to decompose the reaction into a prompt stage, 
in which high momentum Z=l, Z=2 fragments and neutrons are emitted, and a second 
stage involving the decay of the remnant left after this staget3. The analysis presented 
here is based on about 32,000 fully reconstructed MF events for which the total charge of 
the reconstructed Au system was found to be 79±4. Average fragment mass values for a 
given Z were determined and used to reconstruct the MF final state mass A' of the charged 
fragments. The number of free neutrons in the MF final state is used in the determination 
of E* by energy balance. 

The remnant resulting from the prompt stage can be characterized by Z, A, and E*. We 
follow previous practice and report E* on a per nucleon basis. The determination of these 
quantities has been described in detail elsewhereta. Fig. 2 shows the mass distribution of the 
remnants. 
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FIG. 2. Mass distribution of the experimental remnant. 



The most probable mass is A~190. However, the distribution is broad and extends 
down to A ~ 100. The E* distribution is shown in Fig.3a. The distribution peaks at ~2 
MeV/nucleon but extends beyond 16 MeV/nucleon. 

It should be noted that the data were obtained with a minimum bias trigger that elim- 
inated some events with very low E*. These events do not lead to MF. Fig. 4 shows the 
variation of E* and A with the total charged particle multiplicity m. These quantities vary 
in the opposite way with m, with E* increasing with m, and A showing a concomitant 
decrease. At a given m there is a distribution in the values of A and E*. For a given m, 
the width of the distribution in A increases from ~ 1% to ~ 13% while that of the E* 
distribution is ~ 25% over most of the multiplicity rangaB. 
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FIG. 3. Distribution of remnants as a function of (a) total excitation energy and (b) input 
excitation energy to SMM. 



B. The effect of non-thermal expansion 

We have previously shown that some of the excitation energy of the remnant actually 
consists of nonthermal expansion energy, i^^O'Ej. The model input energy Ei is obtained 
by subtracting E^ from the excitation energy E*. Standard SMM does not include radial 
expansion and so only Ei must be used in the input data. Fig.4 shows the dependence 
of E^ and E^ on the multiplicity m. The expansion energy is very small for m^20 but 
then increases sharply, becoming comparable to the input energy for the largest observed 
multiplicity, m ~ 60. The spectrum of E^ values is plotted in Fig. 3b. The distribution 
peaks at slightly less than 2 MeV/nucleon and extends to ~10 MeV/nucleon. 
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FIG. 4. Dependence on total charged particle multiplicity of some average properties of the 

experimental remnants: (a) mass, (b) total excitation energy (E*), nonthermal expansion energy 
(Ex), and input energy (Ei). All energies are in MeV/nucleon. 

The expansion energy was obtained as the difference between the sum of the measured 
charged fragment mean kinetic energies and the translational thermal and Coulomb energies 
of the fragment sli^'Ej. This procedure involves the use of the Albergo double isotope ratio 
thermometerEI to determine the isotopic temperature Ta- This temperature determines the 
translational thermal energy per particle, which is 3/2 T^. As shown in the next section, 
SMM then independently predicts the observed double isotope ratio. This self consistency 
supports the validity of the Ei determination and permits a combined definitive test of the 
isotope ratio thermometer T^ and SMM. 

Additional independent evidence for the presence of expansion energy in the data can 
be seen in a comparison of the mean transverse kinetic energies of fragments with the SMM 
predictionsca. Fig. 5 shows the results for Li-N fragments for five multiplicity bins: 1-9, 
10-19, 20-29, 30-39, and 40-59. Generally good agreement is obtained for low multiplicities, 
as expected, because of the small contribution of the radial expansion energy for small m 
(see Fig. 4). The SMM transverse kinetic energy values decrease with increasing m. 
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FIG. 5. Mean transverse kinetic energies of fragments as a function of nuclear charge. 

multiplicity bins are: 1-9, 10-19, 20-29, 30-39, and 40-59, respectively. 
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In general, the SMM transverse energies are smaller than the experimental values. The 
trend in the SMM values is the result of two factors, both of which lead to a decrease in 
the translational kinetic energy with increasing multiplicity: (1) the increase in the volume 
occupied by the fragmentsE3 and (2) the decrease in the average charge of the fragments. In 
particular, see Fig. 11, which shows the variation of the calculated and experimental charges 
of the largest fragment with m. SMM is in excellent agreement with experiment indicating 
that the reason for the discrepancy shown in Fig. 5 does not lie in the determination of the 
mutual Coulomb energy. Rather, the increasing contribution of non-thermal expansion to the 
experimental energies dominates the change in Coulomb energy. A similar discrepancy has 
been observed when another data set was used to compare SMM with the Berlin statistical 
modela. Here too, the discrepancy was attributed to radial expansion energytj. 

IV. COMPARISON OF EXPERIMENTAL DATA WITH SMM 

A single SMM calculation was performed for each of the 32,000 EOS events. The input 
data consisted of the Z, A, and Ei values of the remnants. The output of each SMM 
calculation gives the Z and A values of each fragment in its asymptotic cold state. The 
distribution of these MF products is used to make the comparisons presented in this section. 
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FIG. 6. Variation of second stage multiplicity with input energy. 



A. Fragment yield and multiplicity comparisons 

Owing to the important role of multiplicity in the EOS experiment, we compare the SMM 
and experimental multiplicities at the outset. This comparison is made as a function of Ei. 
See Fig. 6. Since SMM does not include a prompt first stage, the SMM charged particle 
multiplicity is compared to the experimental second stage multiplicity, m2- This quantit 
is obtained from m by subtracting event-by-event the prompt first stage multiplicity, mil 
The calculated and experimental m2 distributions are in close agreement. Because of this 
agreement, SMM results will be plotted as a function of m in many subsequent comparisons 
with data. Here m will be the sum of SMM 1712 and experimental mi values. 
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FIG. 7. Yield of fragments averaged over all m as a function of fragment charge. 

The fragment yield distributions offer the most direct comparison between the data and 
SMM. Fig. 7 shows the number of Z=3-16 fragments per event averaged over all multiplicities. 
Good agreement is obtained for nearly all the fragments and the overall trend of decreasing 
yield with increasing charge is well reproduced. Both data and SMM show that the yield 
of fluorine is suppressed relative to that of neighboring fragments. SMM indicates that this 
low yield reflects the influence of final state interactions on the primary fragments. 

The individual fragment yields are plotted as a function of m in Figs. 8 and 9. SMM 
generally does an excellent job of reproducing the data except for the lightest fragments, 
where the model predicts too many fragments for large m. Both data and SMM show 
the characteristic rise and fall of the individual fragment yields as a function of m or Ei. 
Both of these distributions were used in the statistical analyses which extract the power law 
and critical scaling behavior. The average value of the total number of intermediate mass 
fragments < IMF > as a function of m is shown in Fig. 10. 
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FIG. 8. Average yield of fragments, for Z = 2, 3, 4 as a function of m. 

The yield of IMF's increases to ~4.4/event at m ~ 48 and decreases for larger m. This 
behavior is reproduced by SMM, although the number of IMF's at the peak is somewhat 
overestimated. This difference is a consequence of the above discrepancy between data and 
SMM for Li and Be fragments at large m. Fig. 11 shows the average value of the charge 
of the largest fragment in the distribution, Zmax- SMM gives the best agreement with the 
data using eo = 16 MeV. 
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FIG. 9. Average yield of fragments, for Z=6, 9, 12 as a function of m. 

We have previouslyOJ obtained from the data the isotopic temperature TqEI on the basis 
of ^H/^H and ^He/'^He double isotopic yield ratios. A value of Ta was also obtained from the 
^Li/^Li and ^He/^He yield ratios. Although the two values of Ta are nearly equal, the TheDT 
values were found to be more robustliS. These temperatures are compared with the SMM 
values in Fig. 12. Excellent agreement in the MF region is observed further confirming the 
self-consistent nature of data and the predictions of SMM. 

It should be noted that the determination of freeze-out temperatures obtained by the 
double isotope yield ratios is subject to correction due to formation of these isotopes in 
secondary decay. This subject has been investigated by a number of workersE3'E2rM. For the 
ThcDT thermometer the correction has been reported ~ 10% below Ei r^ 7 MeV/nucleont3. 
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FIG. 10. Average yield of intermediate mass fragments for Z=3-30 as a function of m. SMM 
results are shown for £„ = 16 and e^ = 8. 
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FIG. 11. Charge of the largest fragment as a function of m. SMM results are shown for eo 
16 and €o = 8. 
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B. The caloric curve using isotopic temperatures 

The asymptotic caloric curve, which is a plot of fragment isotopic temperatures versus 
input energy, was first obtained by Pochodzalla et al. for 600 A MeV Au + Au collisions^. It 
was found that the temperature had the essentially constant value of ~5 MeV for excitation 
energies between 3 and 10 MeV per nucleon, a result that was interpreted as evidence for 
a first-order phase transition. A similar analysis of the EOS data, in which the excitation 
energy included the expansion energy, showed that the temperature increased continuously 
but slowly with E* over the above range, i.e. from ~4 to ~6 MeVll3. Having established 
the presence of non-thermal expansion, we can redetermine our caloric curve and compare 
it with SMM. 
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FIG. 12. Dependence of the ^H/^H to '^He/^He (ThcDt) isotopic ratio temperature on m. 



These SMM and EOS caloric curves are compared in Fig. 13, where both SMM isotopic 
temperatures and experimental isotopic temperatures are plotted versus the experimental 
input energy Ei per nucleon. The SMM and EOS curves are in close agreement and show 
a very sharp increase in the temperature T over the experimental energy range. Included 
for comparison is the caloric curve previously obtained with the inclusion of the expansion 
energy, which shows a much slower variation of T with E*. The determination of a caloric 
curve in the manner proposed in Ref.t3 has already been shown to be problematicHi^. The 
recent reanalysis by the ALADIN group is in close agreement with our datao. 

It must be noted that experimental caloric curves are actually inadequate measures of 
the thermodynamic caloric curve, which involves the breakup temperature rather than the 
isotope ratio temperature. This has been discussed both by Bondorf et al.a and in a recent 
publication from our groupEJ. 
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FIG. 13. The caloric curve using the isotopic temperature Ta = ThcDT- Experimental results 

are shown as a function of both the input energy/nucleon and the total excitation energy/nucleon. 
The SMM results are for the input excitation energy/nucleon. 



C. Critical exponents and related quantities 

The EOS collaboration has analyzed the lA GeV ^^^Au data in terms of the theory of 
critical phenomena, according to which MF is viewed as a continuous phase transition. In 
this section we subject the SMM events to this same analysis and compare the results with 
those obtained from the data. 
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1. The critical point multiplicity rric and the exponent t 



(a) 



10' 



. 1 2.6 't (b) 
• S2.4 - 




10 

FIG. 14. Determination of r and rric from datauil. (a) Xu values from the power law fit to the 
fragment mass yield distribution obtained for different m. Values of r as a function of m. (c) 
Power law fit to data point m = rric, corresponding to the minimum value of xl- T^i^^ dashed line 
is a fit to the open points which exclude A = 4 fragments. The block dot results include the A = 
4 fragments. 

In order to extract the various critical exponents from the data the location of the critical 
point, which can be characterized by the critical multiplicity, mc, must be determined. We 
have used the method presented in referencell3llJ, in which rric, the critical exponents r, 7 
and a, and the scaling function were obtained from the EOS data. We briefly summarize the 



procedure below. The fragment mass yield distribution, A^^i (e), where e 



[rrir 



m)/mc 



IS 



normalized to the mass of the remnant, Aj.em{^)- The normalized fragment distribution can 
be written as^fl 



riAfie) = NAf{e)/Arem{^) = qoAff{z) 



(7) 



where f{z) is the scaling function and the scaling variable z = ^A'f. If we assume that 
scaling is valid for clusters of all sizes, then qo is a function of only r: 



Qo 



1/CO 



11 



(8) 



where ( is the Riemann zeta function and 2 < r < 3l3'l3. At the critical point f{z) = 1 
and a pure power law is obtained for the fragment mass yield distribution. The power law 
behavior is modified by the scaling function away from the critical point. 
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FIG. 15. Determination of r and rric from SMM. See Fig. 14 for details. 

In order to determine m-c and r one performs power law fits to n^ (e) versus Aj for each 
value of m over a particular range of fragments, e.g. Zj = ?> — 16. The best fit, as determined 
by the minimum value of x^; gives the value of m = rric. The results of this analysis for the 
EOS datafij are shown in Fig. 14. The fitting procedure yields a deep minimum in x^ at 
m = rric = 22 ± 1, and thus accurately determines the location of the critical point. The 
power law fit to the data at this multiplicity is also shown. The value of r is 2.19±0.02. 

The above analysis was also performed for the SMM events for Zj = 3 — 16 and the 
results are shown in Fig. 15. The dependence of xl aiid r on m and the quality of the power 
law fit at rric are similar to that exhibited by the data. The results are summarized in Table 
|l| and are in excellent agreement with the data. 

2. Determination of the exponent a 

Once rric is known, it is possible to determine the value of a. We have used the percolation 
procedure where the largest piece Z^ax, A^ax is removed only on the liquid sideE^. 

Since the scaling function must have a single maximuml^, we define z = Zmax = constant 
as the value of the scaling variable where the maximum for this fragment mass occurs. The 
previously given relation between z and a can then be written for this maximum as 






constant 



(9) 



where emax = {ttt-c — fnmax{A.f))/mc, with mmax{A.f) being the multiplicity for which the 
maximum yield of fragments of mass Af is obtained. Eq. (7) leads to a power law, emax oc 
Aj'^ , from which a can be found. 

Applying this analysis to SMM , we first determine the values of mrnaxi^f)- Typical 
results have been shown in Figs. 8 and 9. Fig. 16 shows the power law plot, ln{mmax ~ fnc) 
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versus ln{Af), from which we obtain a = 0.63 ± 0.08. Applying the same analysis to EOS 
data, we obtain a = 0.32 ± 0.05, much smaller than the SMM valued. The analysis in Sec. 
V suggests that the a values in data could be affected by the hot fragment cooling process, 
and that the a value adjusted for cooling would be a = 0.54 ± .11 (See Table II) This value 
is used in constructing the scaling function for the EOS data. 
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FIG. 16. Determination of a from SMM. A linear fit of ln{mmax — i^c) vs ln{Af) for Z=3-15 
fragments gives a = 0.63 it .08. 



3. The scaling function 

Knowing the values of 171^, t, and a for data and SMM, it is possible to evaluate the 
scaling function by rewriting Eq. (5) as 



f{z) = nAf{t)/qoA~. 



(10) 



The unsealed data is shown in Fig. 17. Scaling n^^(e) according to Eq. (8) collapses the 
multifragmentation data from a broad range of fragments yields onto a narrower band for 
both data and SMM, as shown in Figs. 18 and 19. Both scaling functions show a comparable 
scatter of the fragments. By definition, both have a value of unity at the critical point, and 
both have a maximum of comparable magnitude for virtually the same value of the scaling 
variable. 
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point. 
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4- The exponent 7 and moments of the fragment yield distribution 

The exponent 7 has been determined by means of the 7-matching technique as apphed 
to the second moment, M2(e), of the fragment yield distribution^. The second moment is 
defined as 



Mo(e) 



T^ua, 



{^A) 



(11) 



and a similar expression may be used if the fragments are characterized by their nuclear 
charge instead of mass. Again, following the percolation theory procedure, we omit the 
largest fragment from the summation in Eq. (9) only on the "liquid" side of the critical 
point (m < iTLi). All fragments are included on the "gas" side {m > rric). The exponent 
7 is obtained from the power law M2 oc |e|~'^ by searching for various multiplicity regions 
that yield values of 'jgas and '^uquid that agree with each other. The value of iric is again 
determined in the same fit. 
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FIG. 19. The scaling functions for SMM. The intersection of the soUd lines marks the critical 
point 

The results for SMM are shown in Fig. 20. The procedure has also been described in 
Ref.t3 and the results obtained from an analysis of the EOS data are summarized in Table 
|T[ The same analysis can be applied to the SMM events and the results are also summarized 
in Table 0. The experimental and calculated 7 values are consistent within the limits of 
error. 

We have also evaluated several other quantities that have been associated with critical 
behaviorn^'El'Ea. They include the fluctuations in the size of the largest fragment and in 
M2, the magnitude of the peak in the combination of moments 72 = M2Mo/M^, and the 
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determination of r from a plot of InMs versus lnM2. These quantities, too, are in good 
agreement with experiment. 
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FIG. 20. Second moment of the SMM generated fragment mass yield distribution. The largest 
fragment has been removed only from the "liquid" side (m < r7T,c).All points are from SMM. The 
solid points from SMM have been used to determine 7. The solid curve is the fit of the power law 

M2 = C±|e|-T. 

The fluctuations of the size of the largest fragment and of the moments provide an 
independent method to locate the critical multiplicity rric of the MF transition. The first 
moment Mi can also indicate how the nucleons are distributed into light particles, IMFs 
and the largest piece in MF. This distribution can identify the phases present in MF and if 
the coexistence of liquid and gas phases occurs in MF. (See Sec. VI.) 

V. EFFECTS OF COOLING AND OF THE USE OF A SINGLE REMNANT 

The SMM calculations described in the preceding sections give results for deexcited sec- 
ondary fragments formed from a distribution of remnants. These remnants are grouped 
according to multiplicity and we have noted that at a given m there is a distribution of 
remnant charges, masses, and excitation energies. An ideal statistical analysis of multifrag- 
mentation would involve primary fragments formed from the breakup of a unique remnant 
{A, Z) as a function of remnant excitation energy or temperature. In this section we use 
SMM to examine the extent to which the departures from this ideal situation affect the 
values of the critical exponents and related fragment properties. We introduce the terms 
SMM/iot and SMMcow to designate results obtained from SMM for primary fragments and 
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deexcited fragments, respectively. 



A. Hot and cold SMM fragment yields 

The most direct view of the effect of deexcitation is provided by a comparison of the 
SMMhot and SMMcojd fragment yield distributions. For simplicity, we evaluate these distri- 
butions for the decay of a single remnant as a function of Ei. We have chosen the remnant 
formed in the Au on C interaction corresponding to the critical multiplicity in the data, 
A = 160, Z = 64, and evaluated its breakup for a range of E^ values. At a given m the 
total rms width of the remnant distribution is ~ 5-7 %. In SMM the average Z/A ratio of 
the hot "pre- fragments" is the same as the remnant Z/A ratio, i.e., no particles or energy 
leaves the remnant system during the ~ 100 fm/c time frame for MF. 
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FIG. 21. Second stage multiplicity m2 for SMMcoid and SMM/jo^. 



The increase in the charged particle multiplicity 1712 from SMMhot to SMMcoid is shown 
in Fig. 21. This figure shows that cooling is important. Typical results for fragments with 
Z > 3 are shown in Figs. 22 and 23. The fragment multiplicities for SMM^o^ and SMMcoid 
are in close agreement up to E^ ~ 7 MeV/nucleon. However, at higher excitation energies 
the primary fragment yields are significantly larger than those of the corresponding cold 
fragments. At these high Ei, fragments undergo substantial secondary decay, which accord- 
ing to SMM, occurs by evaporation for fragments with A > 16. For A < 16, secondary decay 
occurs by Fermi breakup. These mechanisms tend to form the lightest particles and thereby 
reduce the yields of fragments with Z >3. We emphasize the fact that the net effect of these 
secondary processes on the fragment distributions is rather small below 7 MeV/nucleon, the 
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energy range where most of the MF occurs. The critical point corresponds to a much lower 
excitation energy, Ei ~ 4.3 MeV/nucleonlla, and therefore quantities determined at these 
low energies should be nearly unaffected. 
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FIG. 22. Multiplicity of Z=3, 4, 5 fragments as a function of input energy for SMM^ot and 
SMMcoid in MF of A = 160, Z = 64. 

Fig. 24 shows the results for Z=l and Z=2 particles. Contrary to the results shown 
in Figs. 22 and 23, most of these particles are formed in the cooling step. Averaging over 
the excitation energy range important for MF, 3-7 MeV/nucleon, we find that the lightest 
hot fragments have a higher internal energy per nucleon than the heavier hot fragments, ~4 
MeV/nucleon versus ~3 MeV/nucleon. Consequently the main source of these particles is 
the Fermi breakup of the lighter fragments. 
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FIG. 23. Multiplicity of Z=8, 10, 12 fragments as a function of input energy for SMM/^ot and 
SMMcoZd in MF of A = 160, Z = 64. 

B. Effect of cooling on critical exponents and related quantities 

The analysis described in Section IVC has been applied to the products of SMMhot, 
using the experimental remnants as input. The resulting values of the various exponents are 
tabulated in Table |ITl| . The calculated value of r for hot fragments is ~ 6% smaller than 
that for cold fragments while that of 7 is ~ 11% smaller. In contrast, the determination of 
the exponent a differs significantly for the two distributions, with the SMM/io^ value being 
~40% larger than that from SMMco^d. The large difference in a values is understandable 
given that a is determined in part by the multiplicities for which the various fragments 
attain their highest yields. As shown in Figs. 8 and 9, the yields of the lighter fragments 
peak at large values of m, corresponding to Ei values of 7 MeV/nucleon and higher. It is 
precisely at these energies that secondary decay reduces the yield of these fragments and 
thereby leads to a smaller a. 
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FIG. 24. Multiplicity of Z=l and 2 particles as a function of input energy for SMM/jo^ (solid 
circles) and SMMcoid (open circles) in MF of A = 160, Z = 64. Note that most of these particles 
are produced in the cooling process. 

In a recent publication Mastinu et al.c3 used the Copenhagen model to examine the effect 
of secondary fragment decay in the MF of ^^^Au. In particular, they compared the shape 
of the second moment, M2, of the fragment charge distribution for SMMhot and SMMcoid- 
They found that the shape of the M2 distribution in the vicinity of its peak is unaffected 
by secondary decay when plotted against E* whereas there is a change in shape when M2 
is plotted against m. 

On this basis they concluded that secondary decay affects the values of exponents such 
as 7, which is obtained from M2, when charged particle multiplicity is used as the control 
parameter. Unfortunately, they did not calculate the magnitude of this effect. We have 
already evaluated the effect on 7 of secondary decay and, as shown in Table |T|, find that it 
is a relatively small 10% correction. We note in passing that the M2 plot shown by Mastinu 
et al.c3 is not appropriate for the extraction of 7 because the largest fragment has been 
removed from both the liquid and gas sides instead of only on the liquid side, as was done 
in Fig. 20. Furthermore, they did not turn off the fission channel in their SMM calculation 
and the contribution of this mechanism disproportionally affects the analysis. 

C. Critical exponents from the multifragmentation of a single remnant for SMM^ot 

The analysis described in Section IVC is applied to the breakup of a single remnant, 
using the input energy as the control parameter. We use the critical remnant A = 160, 
Z = 64 as the system size for SMM. Fig. 25 shows the SMM/^ot results for the determination 
of the critical point by means of the one-parameter power law fit to the fragment mass yield 
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distribution. 

Again, there is a deep minimum in xl which determines the location of the critical point. 
This should be compared with the same results for SMMcoid using the experimental remnant 
distribution and m as the control parameter (Fig. 15). 




10 Af 

FIG. 25. Determination of r and E^ in the MF of A = 160, Z = 64 for SMMhot- The panels 
have the same meaning as in Fig. 14 except that Ei replaces m as the control parameter. 

The value of r obtained from the fit in Fig. 25 is given in Table JTT[ Also included 
are the values of the other exponents. A comparison with the corresponding exponents 
obtained for SMMcoid for the experimental remnant distribution indicates that the difference 
in r and 7 values is less than 10%. Thus the combined effect of secondary decay, the 
remnant distribution for a given m, and the use of multiplicity as the control parameter 
have only a minor effect on the values of the r and 7 exponents. However, the exponent a 
for SMMhot is much larger than the corresponding value of SMMcow just as was the case for 
the experimental remnant distribution. Again, this difference reflects the effect of cooling. 
We note here that the use of a single remnant and excitation energy as the control parameter, 
which is most readily seen in a comparison of the two SMMhot columns, has essentially no 
effect on the critical exponents. 
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FIG. 26. Unsealed yield of fragments/event from SMM/^ot 
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for A 



30 as a function of 



Fig. 26 shows the unsealed SMM/io^ fragments for A = 5 to 30 for the single remnant 
system. Fig. 27 shows the scaling function obtained from these data and the resulting scaling 
collapse of the data into a very narrow band. This scaling function should be compared with 
the scaling function from SMMcoZd (Fig- 19). The more extensive scatter in Fig. 19 can be 
attributed in comparable measure to the effects of fragment cooling and to the presence of 
the remnant distribution, using multiplicity as the control parameter. This indicates that 
the scatter seen in the experimental data at the ends of the scaling curves, Fig. 18, may 
also be due to these effects. 
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VI. NATURE OF THE PHASE TRANSITION IN SMM 

We have shown that SMM can reproduce the various features of the EOS data, including 
the values of most of the critical exponents, critical multiplicity, and critical scaling. If 
SMM did not predict the occurrence of a continuous phase transition then the exponents 
would merely constitute a particular parameterization of the data and would have no ulterior 
significance. It is therefore important to determine if SMM predicts such a thermal phase 
transition. 

We follow the approach of Grosso and HiilleiO, who show that the nature of the transi- 
tion can be determined by means of the microcanonical equation of state, using a plot of the 
reciprocal partition temperature j3p = l/Tp versus Ei. We assume Tp is the best estimate 
of the average value of the fluctuating event-by-event SMM microcanonical temperature. 
A first order transition is identified by backbending in the jSp vs Ei curve and hence has 

,d(3p 



a region of negative system specific heat, C„ 



^P' dEi 



A continuous phase transition 



will not exhibit backbending and will have a positive system specific heat that peaks at the 
critical energy. 
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FIG. 28. Dependence of fip on the input energy for A=70, 100 and 130. The lines are drawn 
for equal area and show the "Maxweh Construction" . 



Recall that for infinite neutral nuclear matter a critical phase transition is expected at a 
temperature ~15-18 MeVc^. At lower temperatures, infinite neutral nuclear matter would 
exhibit a first order liquid-gas phase transition with possible coexistence of constant density 
phases. SMM, however, describes the transition of a highly excited finite charged nuclear 
drop into a number of hot intermediate mass fragments followed by a de-excitation process. 
The cold state liquid and gas consists of smaller droplets and nucleons emitted in the cooling 
process. In the hot stage, the so-called pre-fragments are produced in thermodynamic equi- 
librium where the volume, surface. Coulomb, symmetry and translational terms determine 
the fragment yieldg23. 

To investigate the interplay of these energy terms we study SMM/ioi as a function of the 
remnant mass and chargeO. Fig.28 shows (3.p versus Ei curves for A = 70, Z = 30; A = 100, 
Z = 40; and A = 130, Z = 53. A characteristic backbending is observed for these systems. 
Following the microcanonical prescription, we make the "Maxwell construct ion" l3, which 
determines the average transition temperature and the magnitude of the transition energy 
or effective latent heat. For A=100 we obtain a transition energy of 3.2 MeV/nucleon as 
given by the width of the backbend region. Figure 29 shows that the system specific heat 
for A = 100 is negative indicating the possibility of a first order phase transition. We have 
computed the effect of the variable volume on the system specific heat and found that it is 
a very small fraction of the total. The translational energy contribution is also very small. 
Thus Cn describes the change in the internal energy of the MF system. The transition 
energy plays the role of an effective latent heat for MF. SMM/jo^ indicates that the effective 
latent heat is reduced for heavier remnants. For the A = 70, 100, and 130 systems, there is 
a progressive decrease in the transition energy as shown in Fig. 30. A linear extrapolation 
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suggests that the transition energy as an effective latent heat could vanish for A > 170. 
We can gain further insight about these transitions by studying the distributions of light 
particles and fragments prior to the de-excitation or cooling process. For A = 100 the first 
moment of the hot mass distribution as a function of the input energy is shown in Fig. 31. 
For Ei = 8 MeV/nucleon only ~ 5 out of 100 nucleons are not part of the A>4 droplet 
distribution. Thus hot SMM describes a transition where a large droplet is multifragmented 
into smaller droplets. The subsequent deexcitation processes are responsible for most of the 
light particles with A <5 seen in the asymptotic mass distributions. 
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FIG. 29. System specific heat per particle for A=100. The solid curve is to guide the eye. 

In this SMM model backbending can be a relatively large effect. The extra surface 
energy required for the MF transition cools the system and the temperature falls. For a 
larger system the fractional change in the surface to volume ratio due to MF is reduced 
and the backbending is smaller. In SMM the total Coulomb energy is reduced by the factor 
1 — 1/(1 + kY^^ ~ 0.3 due to expansion prior to clusterization. Since the total energy is 
measured with respect to the ground state of the remnant system A, Z at normal density, 
this introduces a Z dependence into the Pp versus E^ equation of state. 

Fig. 30b shows the average transition temperature, which decreases as the remnant mass 
and charge increases. This decrease clearly reflects the increase in Coulomb energy for the 



34 



heavier reninantgI3. In contrast, for a liquid-gas phase transition in finite neutral matter 
we would expect the transition temperature to increase with the system size, because the 
binding energy per nucleon increases. For a larger charged nuclear remnant, the transition 
energy will vanish if enough Coulomb energy is converted to heat energy in the expansion 
of the remnant prior to the MF transition. 
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FIG. 30. a) MF transition energy, the effective latent heat, for A = 70, A 
b) Average, i.e., Maxwell construction, transition temperature for A = 70, A 
A = 160. 



100, and A = 130. 
100, A = 130, and 



For the remnant mass A = 160 the reciprocal partition temperature j3p is plotted versus 
Ei in Fig. 32a and the fluctuating event-by-event reciprocal microcanonical temperatures 
f3p, are shown in Fig. 32b. The area under the j3p vs Ei curve is the system entropy. In 
Fig. 33a, we expand the Pp vs Ei plot. The system specific heat per particle is computed 
by taking differences between adjacent points. C„ is shown in Fig. 33b. The specific heat is 
positive and peaks near the critical energy, indicating a continuous phase transition. Hiiller 
characterizes such behavior as a continuous phase transition with an anomalous specific 
heat!§ 
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FIG. 31. First moment of the SMMhot mass distribution for A = 100. 

A method to determine the order of the phase transition based on canonical kinetic 
energy fluctuations has recently been proposed&ilj. To test the conclusions based on the 
SMM/ioi microcanonical equation of state, we have done this type of analysis on the Au on 
C data. For a single remnant system size A, Z this analysis compares the canonical specific 
heat Ci with the thermally scaled kinetic energy fluctuations {E'^— < E >^)/T^. 

Following the prescription used on the 35 MeV/nucleon Au on Au data we find that the 
energy fluctuations from the Au on C data are dramatically reduced when a AA/A < 0.03 
remnant mass cut is made, which assures that we have a single system size, the necessary 
constraint for the fluctuation analysis. Thus for the Au on C data, we flnd that Ci > 
{E"^— < E >^)/T^ in the whole MF region. The value of Ci still remains positive even for 
a larger remnant mass cut e. g. AA/A ~ 10%. This result argues for a continuous phase 
transition. A report of this work is in preparation. 

Very recently, (since the submission of the present article) theoretical analyses of the 
nature of MF phase transition and its effect on the caloric curve have been publishedE3'E3, 
indicating the occurrence of a flrst order phase transition. Both of these calculations neglect 
the Coulomb interaction, which has been shown to play an important role in MFII3. It is 
obvious that there will be a flrst order phase transition in absence of Coulomb energy. 
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FIG. 32. a). Dependence of Pp on the input energy Ei for A = 160. b). Event to event 

fluctuations of /? versus the input energy for A = 160. The dot scatter points identify compound 

nucleus events, the box scatter points identify MF events. The soHd block dots are the average /3 

values for MF events. 

It has been suggested on the basis of a caloric curve using the SMM partition temperature 
versus an experimentally determined input energy, that MF of gold is a first order liquid-gas 
phase transitionO. This curve, computed for comparison with the ALADIN data showed 
that the temperature first increases with excitation energy, then remains nearly constant at 
5-6 MeV as Ei increases between 3 and 10 MeV/nucleon, and then again increases. The 
first regime has been interpreted as the liquid phase, i.e. the compound nucleus, the second 
as the MF coexistence phase, and the third as the gas phase, consisting of a mixture of 
nucleons and a few of the lightest fragments^. This viewpoint can be probed by examining 
the first moment of the hot fragment yield distribution for a heavy remnant. Fig. 34a shows 
a plot of the dependence on E^ of the first moment of the hot fragment yield distribution 
for the experimental remnant distribution. 
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FIG. 33. a). Expanded (3p vs Ei curve, b). System specific heat per particle for A=160. 

The figure shows that even for Ei ^ 8 MeV/nucleon, well above the MF transition, only 
about 8% of the remnant mass ends up in particles with A < 4 and that the other 92% 
are in intermediate mass fragments with A >4. This argues against the coexistence of the 
constant density liquid and gas phases in SMM. For the experimental remnant distribution, 
we compare the first moment Mi, for SMM/jo^ Fig. 34a with Mi for SMMcoid as shown 
in Fig. 34b. The cooling of the hot fragments which carry the multifragmentation signal 
produces a large number of final state nucleons and light composites in agreement with the 
data. The IMFs survive the cooling process and identify the MF transition. This is due in 
part to the low ~ 4.3 MeV/nucleon critical energy for MF. For an A = 70, Z = 30 remnant 
, the predicted center of the backbend region energy is ~ 7.5 MeV/nucleon. Here, the hot 
fragments will have to emit a significantly higher fraction of their fragment mass in the 
cooling process and the MF signal could be severely attenuated. Thus the Coulomb energy, 
which lowers the MF threshold and reduces the effective latent heat, plays a central role in 
MF. Here the vanishing transition energy reflects the changes in the surface, volume, and 
Coulomb energies associated with MF rather than the traditional latent heat in a liquid-gas 
phase transition, which primarily involves the transfer of individual nucleons from the liquid 
to the gas, where both phases are at constant density puquid and Pgas- 
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FIG. 34. a) First moment of the mass distribution from SMM^oid for the experimental remnant 
distribution, b) First moment of the mass distribution from SMM^ot- The coohng mechanism is 
nucleon and A < 4 composite particle emission. 



VII. SUMMARY AND CONCLUSIONS 

We have compared the EOS inultifragmentation results for lA GeV Au on C with the 
SMM model where the volume parameterization is determined by experiment. We have 
found that the standard SMM parameters produced excellent agreement. The input to 
SMM consisted of the Z, A, and Ei of the individual remnants produced in some 32,000 fully 
reconstructed events. Thus, in contrast to earlier experiments which have been compared 
with SMM, the A, Z, Ei event-by-event input data were obtained from experiment rather 
than by using SMM to constrain average input conditions, or by the use of a theoretical first 
stage transport model. 

We find that SMM is in very good agreement with the observed fragment charge yield 
distribution, fragment multiplicity distributions, total charged particle multiplicity, isotope 
ratio temperatures and caloric curve. This agreement is obtained when the expansion energy 
is subtracted from the experimentally determined excitation energy of the remnant. The 
agreement provides further confirmation that multifragmentation can be described as an 
equilibrium process. 

SMM also predicts a power law for the fragment mass yield distribution at essentially 
the same multiplicity or excitation energy as is observed in the data. Other features charac- 
teristic of a continuous phase transition, such as the scaling function, are also reproduced. 
The critical exponents r, 7, and a were obtained from SMM for comparison with the exper- 
imental values. The agreement is excellent for r, which is based on results obtained just at 
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the critical point, fair for 7 which is based on results obtained above and below the critical 
point, and poor for a, which is based on results well above the critical point. 

The effect of cooling of the primary hot fragments produced by SMM has been evaluated. 
The primary hot fragment yield distributions for Z >?> are only minimally affected by the 
cooling process below Ei =7 MeV/nucleon. The SMM a values are substantially affected by 
cooling, while r and 7 are unaffected by cooling. The scaling function obtained for SMM/ioi 
is fully collapsed into a very narrow band. In contrast to fragments with Z>3, the yield of 
particles with Z<2 is substantially increased by cooling, primarily as a result of the Fermi 
breakup of light fragments. 

We have also performed SMM calculations in which the experimental remnant distribu- 
tion has been replaced by a single average remnant. The results are virtually unchanged 
indicating that the mixing of different remnants in the experimental data does not affect 
the results. We have used both multiplicity and excitation energy as the control parameter 
in the determination of critical exponents and related quantities and find no difference. 

The nature of the phase transition predicted by SMM has been examined using the 
microcanonical equation of state. For lighter remnants we find evidence of backbending 
in the caloric curve and a negative system specific heat, which are the signatures of a 
first-order phase transition. We estimate the transition energy by means of a Maxwell 
construction. The transition energy decreases with increasing remnant mass and charge 
and may extrapolate to zero just above A = 160, which might suggest a continuous phase 
transition in the breakup of this remnant. Here the positive system specific heat peaks at 
the same Ei value for which both data and SMM fragment mass yield distributions obey a 
power law. For both the first order and continuous phase transition cases, SMMhot indicates 
that the multifragmentation final state consists of droplets with A > 4 and that particles 
having A < 4 account for only ~ 8% of the mass in the MF region. 

The SMM results agree with the theoretical expectations for a small (~ 150 constituents) 
isolated systemB^'E^IIa. The most probable equilibrium state of a highly excited small iso- 
lated system with short range interactions is an inhomogeneous state. The addition of the 
long range Coulomb force lowers the MF transition temperature significantly for heavier 
remnants and can influence the order of the thermal phase transition. Both the micro- 
canonical equation of state and thermally scaled kinetic energy fluctuation arguments favor 
a continuous phase transition. 
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TABLES 
TABLE I. Parameters in the SMM free energy expression. 



Fixed Parameters 



Wq Volume binding energy of cold nuclear matter 16 MeV 

(3o Surface tension of the cold nucleus 18 MeV 

Tf, Neutral matter critical temperature 18 MeV 

^sym Symmetry energy coefficient 25 MeV 

d Crack width parameter 1.4 fm 

K Coulomb reduction parameter 2 



Single adjustable Parameter 


eo 




Inverse density parameter 




16 MeV 




TABLE II. 


Critical exponents from data and SMMcoid 


(experimental remnants). 




Parameter 


Data 


SMM 





rUc 22±1 20±2 

r 2.19±0.02" 2.17±0.02 

7 1.4±0.3" 1.02±0.23 

a 0.32±0.05^ 0.63±0.08'= 

0.54±0.11^ 

a. The results differ slightly from those given in Ref. [10] because a larger data set was 

used in the present analysis [58]. 

b. The first value is obtained by including the largest fragment on the "gas" side of mc. The 

second value was obtained by adjusting a for the effect of cooling. The experimental 
value for the scaling function in Fig. 18a was obtained using the second value of a. 

c. This result is obtained by including the largest fragment on the "gas" side of rric. 
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TABLE III. Critical exponents from SMM/,ot and SMM, 



cold- 



Parameter 


Experimental 


remnant 


Sin 


Lgle remnant 




distribution 


A = 


: 160, 


Z : 


= 64 




SMM^oid 


SMMhot 


SMMcoM 






SMMf^ot 


mc 


12±2 


5±2 


12±2 






5±2 


T 


2.17±0.02 


2.05±0.02 


2.07±0.01 






2.03±0.01 


7 


1.02±0.23 


0.91±0.20 


1.10±0.20 






0.97±0 .15 


a 


0.62±0.08 


1.04±0.11 


0.69±0.02 






1.00±0.07 
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